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ABSTRACT 

The gamma-ray burst (GRB) 060218 has ~ 10 5 times lower luminosity than typical long GRBs, and is 
associated with a supernova (SN). The radio afterglow displays no jet break, so that this burst might arise 
from a mildly-relativistic spherical outflow produced by the SN shock sweeping the stellar surface. Since this 
model is energetically difficult, we propose that the radio afterglow is produced by a non-relativistic phase of 
an initially collimated outflow (jet). Our jet model is supported by the detection of optical linear polarization 
in the SN component. We also show analytically that the jet can penetrate a progenitor star. We analyzed the 
observational data of the prompt emission of this burst and obtained a smooth power-law light curve which 
might last longer than 10 6 s. This behavior contrasts with the long intermittent activities with the X-ray flares 
of typical GRBs, implying that the central engine of this burst is different from those of typical GRBs. This 
argument is consistent with the analysis of the SN component of this burst, which suggests that the progenitor 
star was less massive and collapsed to a neutron star instead of a black hole. The collimation-corrected event 
rate of such low-luminosity GRBs is estimated to be ~ 10 times higher than that of typical long GRBs, and 
they might form a different GRB population: low-luminosity GRBs are produced by mildly-relativistic jets 
from neutron stars at the collapses of massive stars, while typical long GRBs by highly-relativistic jets from 
black holes. We suggest that the central engine of GRB 060218 is a pulsar (or a magnetar) with the initial 
rotation period Po ~ 10 ms and the magnetic field B ~ 10 16 G. A giant flare from the magnetar might be 
observed in future. 

Subject headings: gamma rays: bursts — gamma rays: theory — supernovae: general 



1. INTRODUCTION 

The gamma-ray burst (GRB) 060218 is the second near- 



est event (z = 0.033) and it is spectroscopically 
ated with the supernova (SN) 2006aj (Modiaz et al. 


associ- 
[2006] 


Sollerman et al.1 120061: iPian et al. 


2006; Mazzali et al. 


l2006t 


Mirabal et al.l 120061: ICobb et al l 


2006; iFerrero et al.1 


2006) 



Such a long GRB/SN association is common and this event 
further supports the established picture that all long GRBs 
are re lated to the deaths of massive stars (for recent re- 
views. IWooslev & Blooml l2006l lMeszarosll2006l |plranll200l 
Zhang & Meszaros 2004). In the most popular model of long 
GRBs, the so-called collapsar model, the core of a massive 
star collapses to a black hole or a neutron star, which drives a 
highly relativistic jet, breaking out the star and making a GRB 
dWooslevlll993HMacFadven & Wooslevl[l999l The relativis- 
tic speed of the outflow is require d for the non-thermal prompt 
emission dLithwick & Sarill2001l and references therein). The 
collimation of the outflow is strongly suggested by a break in 
the afterglow, s ince it is produced by the sideways expansion 
of the jet (e.g., iRhoadsll 19991 ISari et al.lfT999t lHarrison et ail 
fl999h . 

However, the prompt and aftergl ow emission o f GRB 
060218 have many peculiarities ([C ampana e t al.l 



Soderberg et all 120061: iGhisellini et all l2006bt iButleri 
Liang et al. 20063): 



2006; 



20061 



(1) The duration of the prompt non-thermal emission in the 
high energy band (15- 150 keV) is St 



10 3 s, and thus 
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this event is one of the longest bursts. 

(2) The isotropic-equivalent luminosity of the prompt 
emission is extremely low ~ 10 47 erg s -1 , which is 
about 10 times lower than those of typical cosmolog- 
ical GRBs. The isotropic-equivalent energy, extrap- 
olated to the 1 - 10 4 keV band in the central engine 
frame, is ~ 6 x 10 49 erg. 

(3) The spectrum of the prompt emission is quite soft com- 
pared to those of typical bright GRBs, and the aver- 
aged spectral peak energy in the central engine frame is 
~ 4.9 keV. (Note that this event obeys the Amati corre- 
lation, despite the other peculiarities. See §|2]formore 
details.) 

(4) Thermal components are detected in the X-ray band at 
t < 10 4 s and in the UV/optical band at 10 4 < t < 10 5 s, 
while other GRBs do not exhibit such a clear thermal 
component. 

(5) The X-ray, UV/optical, and radio afterglows show chro- 
matic features. The X-ray afterglow (at t > 10 4 s) de- 
cays with a standard temporal slope, but has a spec- 
trum much steeper than those of typical GRB X-ray af- 
terglows. The UV/optical afterglow is quite dim and 
dominated by the thermal component and the SN com- 
ponent. Only the radio afterglow seems rather typ- 
ical and explainabl e in the standard external shock 
synchrotron model dSari et al. 1998; R ees & Meszarosl 
119921; iPaczvhski & Rhoadsll 19931) . but does not show a 
jet break until t ~ 22 days. 

Although such a low-luminosity event looks rare, the in- 
trinsic event rate could be very high ~ 10 2 Gpc~ 3 yr" 1 
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compared with the local rate of typical long GRBs deduced 
from the BATSE data R hG ~ 1 Gpc" 3 yr" 1 dSoderberg et al I 
2006tlPianetal.ll2006t iCobb et al]|2006t iLiang et al.l l2006b: 
Guetta et ah! 120051) . For this reason, it has been actively 
debated whether low-luminosity GRBs form a new GRB 
population and whether they have intrinsically different out- 
flow mechanisms and emission mechanism s (see Stane k etal] 



2006; Ghisellini et al. 200 6a]; iKaneko et a l. 2006; Amati et al.1 
2006t IStratta et al.ll2006t iDai et al.l l2006t IWang et al.ll2006l) . 
It has been also argued that the high energy neutrino back- 
ground from the low-luminosity GRBs could b e comparable 
with or larger than that fro m typical long GRBs (Mura se~et al.l 
120061 iGupta & Zhangil2006l) . 

For GRB 060218, it has been widely suggested that the out- 
flow is spherical since the radio afterg low has no jet break 
(ISoderberg et al.ll2006t iFan et alj|2006h . If true, the collapsar 
model cannot be applied to this event because the outflow be- 
comes non-relativistic by loading all the matter of a progenitor 
star. The relativistic spherical outflow might be produced by 
the outermost parts of the stellar envelope that the SN shock 
accelerates when propagati ng through the steep density gradi- 
ent near the stellar sur face dColgate)ll974HMatzner & McKed 
1999; iTan et aO200Tl) . However. iTan et al.l (120011) has shown 
that the energy ~ 10 48 erg is transferred to mildly relativis- 
tic material when the kinetic energy of the SN is £sn ~ 
10 52 erg. For GRB 06 0218, E SN is estimated as ~ 2 x 10 51 erg 
(iMazzali et al.l |2006), so that it is quite unlikely that the 
prompt non-thermal emission with E 7 j S o — 6x 10 49 erg is pro- 
duced by this type of the outflow (see also Matzer 2003b. O 
(2006) has also shown that the energy of the thermal com- 
ponents > 10 49 erg is too large to be explained by the shock 
breakout in the underlying SN. 

In this paper, we show that GRB 060218 can be produced 
by the standard collapsar jet model. We show that the avail- 
able radio data may be interpreted as a non-relativistic phase 
of the exter nal shock after the jet break within the standard 
model (e.g.. lFrail et al]l2000t ITivio & Waxmanll2000l) . We ar- 
gue that the outflow with an initial opening angle 9q ~ 0.3 
and Lorentz factor To ~ 5 can produce the synchrotron radia- 
tion which explains the radio afterglow and is compatible with 
the UV/optical and X-ray afterglow (§ |3). We also examine 
whether such a wide and weak jet can penetrate a progeni- 
tor star b y extending ana lytical considerations of the collapsar 
model by lMatzel (120031) (§l4l 

Remarkably, the detection of optical linear polarization in 
the SN component of this event has been recently reported 
dGorosabel et al.ll2006l) . This observation strongly supports 
our arguments that GRB 060218 arises from a jet. 

Within the jet scenario, there is a possibility that 
low-luminosity GRBs arise as typical cosmo l ogical 
GRB jets viewed off-a x is (e.g., lYamazaki et al.1 120031: 



Rami rez-Ruiz et all 120051; iGranot et al 



2005) 
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The relativistic jet emits 7-rays into the expansion 
direction through the beaming effect. Thus, if the jet is 
viewed off-axis, the 7-ray flux is strongly suppressed. This 
scenario leads to similar event rates of typical GRBs and 
low-luminosit y GRBs, which seems inconsistent with the 
observations dCobb et al.l 120061) . Aside from this statistical 
argument, we also examine the off-axis scenario for GRB 
060218 and conclude that this scenario is unlikely for this 
event because an unrealistically high 7-ray efficiency is 
required. 

In addition, we analyzed the data of the Burst Alert Tele- 



scope (BAT; iBarthelmv et al.l l2005) and the X- Ray Telescope 
(XRT; iBurrows et al.ll2005al) of Swift satellite dGehrels et al.l 
2004- We found that the non-thermal component of the 
prompt emission of GRB 060218 may be fitted by the Band 
function with the spectral parameters similar to those of typ- 
ical GRBs. Furthermore, we found that the light curve of 
the non-thermal component evolves smoothly as a power-law 
function of time, and it might last so long duration as to con- 
nect to the anomalous X-ray afterglow detected up to ~ 10 6 s. 
This behavior contrasts with the long intermittent activities 
with the X-ray flares of typical GRBs, which might be a sign 
of the difference of the central engines. This is also suggested 
by the analysis of the SN component of this event. The SN 
spectrum has less broad lines than those of other GRB-SNe 
and lacks oxygen lines, and the light curve evolves somewhat 
faster. All these facts indicate that the total kinetic energy 
£sn and the ejected mass M e \ ar e both less than those of other 
GRB-SNe. IMazzali etail d2006l) performed a detailed mod- 
eling of the spectra and light curve of the SN component and 
obtained £sn — 2 x 10 51 erg and M e j ~ 2M Q . Then they have 
suggested that the progenitor star of this burst was less mas- 
sive than those of typical long GRBs and collapsed to a neu- 
tron star instead of a black hole. 

Our goal is showing that the low-luminosity GRB 060218 
has a prompt non-thermal emission with a typical Band spec- 
trum and may originate from a standard collapsar jet possibly 
driven by a neutron star. We will not discuss the emission 
mechanism of the prompt non-thermal emission and the ther- 
mal emissions. This paper is organized as follows. In § [2] 
we show the results of our analysis of the BAT and XRT data 
and suggest several implications for the nature of the prompt 
non-thermal emission. In §|3] we display a jet model of this 
event, and discuss whether the jet can make a hole in the star 
in §[4] A summary and discussion are given in §|5] Since GRB 
060218 is very close (z = 0.033), we will neglect the cosmo- 
logical effect, i.e., we set z = 0, for simplicity throughout this 
paper. 

2. THE PROMPT NON-THERMAL EMISSION 

2.1. BAT and XRT data analysis 

We performed the analysis of the BAT and XRT data using 
the standard Swift software package (HEAsoft 6.0.4) and the 
CALDB 2005-11-28. The detector plane histogram (DPH) 
data were used to extract the spectra of BAT. Before produc- 
ing the spectra, we applied baterebin to rebin the DPH 
data using the most accurate non-linear energy correction. 
The detector map for disabling noisy detectors in the analysis 
was created by bathotpix. The mask-weighting map was 
created by batmaskwt img using the optical afterglow posi- 
tion. With including the detector map and the mask- weighting 
map to batbinevt, 13 BAT spectra were extracted from 
the DPH file which contains the data just after the space- 
craft slew. The BAT response matrices were generated by 
batdrmgen. The systematic error vectors were applied to 
the BAT spectra using batphasyserr. For the XRT data, 
we obtained the cleaned event file in the window timing (WT) 
mode from the Swift HEASARC Archive. The foreground 
was excluded by the box region of 1.2' x 0.6'. The back- 
ground region was extracted by the same region size, but out- 
side of the foreground source region. 13 XRT foreground and 
background spectra using the same time intervals of the BAT 
spectra were extracted. The auxiliary response file was gener- 
ated by xrtmkar f . The XRT spectra were binned to contain 
a minimum of 20 photons for each spectral bin. 
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The XRT and BAT spectral data were analyzed jointly with 
XSPEC 11.3.2. The energy ranges used in the analysis were 
0.5-10 keV and 14.0-150keV for XRT and BAT, respectively. 
We multiplied the constant factor to the spectral model to 
take into account the calibration uncertainty in the response 
matric es of the instruments. As reported by the several au- 
thors dKaneko et alj |2006t ICampana et aTll2006t lButlerj|2006t 
iLiang et alj|2006d) . the blackbody component was necessary 
to obtain a good fit for the low energy part of the XRT spec- 
trum. Thus, we performed the fitting with an additional black- 
body (BB) component to an absorbed Band function and an 
absorbed power-law times exponential cutoff (CPL) model. 
For the purpose of this paper, we are interested in the spec- 
tral parameters of the non-thermal component. To obtain the 
well constrained spectral parameters of the non-thermal com- 
ponent, we fixed the parameters of the absorption (A^h) and 
the BB component to their best fit values of each fit, and then 
calculate the uncertainties of the spectral parameters of the 
non-thermal component. The results are shown in TableQ]and 
|2] The most of the spectra were well fitted with an absorbed 
CPL model with an add itional BB component as reported by 
iKaneko et all (120061) and lCampana et al.ld2006l) . However, we 
found a significant improvement in a fit with the Band func- 
tion for the last four time intervals from 1550 s to 2732 s after 
the BAT trigger. The differences in \ 2 between the CPL and 
Band fit were 9.6, 1 1.7, 17.0, and 13.8 in 1 degree of freedom 
for these spectra. Based on this result, we decided to adopt 
the spectral parameters derived by the Band function for all 
time intervals. The temperature of the BB component varies 
slightly between 0.12 and 0.29 keV. Thi s behavior is consis- 
tent with the result of the other an alyses dKaneko et al.ll2006t 
ICampana et al.ll2006tlButIeril2006l) . 

Figure Q] shows the result of the temporal variations of the 
low-energy and high-energy photon index, as and (3b- Fig- 
ure |2] shows the light curve in the 15- 150 keV band and the 
spectral peak energy E p as a function of time. At the flux de- 
cay phase (f > 500 s), the light curve and spectral peak energy 
are well described by power-law functions of time, F oc t 
and E p oc f" L6 (dotted lines). 

Investigating t he consistency of the BAT light curves pre- 
sented by us and lCampana et a l. (2006), we found systematic 
differences especially for t > 1000 s. Although both of the 
BAT light curves are bas ed on the joint spectr al analysis of 
the BAT and XRT data, ICampana etail (120061) use the CPL 
model while we use the Band function. We no ticed that the 
fluxes of the light curve of ICampana et al.l(l2006l) at t > 1000 s 
are in the level of ~ 10 mCrab. BAT can detect the Crab neb- 
ula at 6a level in 1 s exposure. If the BAT sensitivity can be 
scaled as a square root of the exposure time ( Mar kwardt et al.l 
l2005h . we need the exposure time of ~ 10 4 s to detect the 
10 mCrab source assuming the Crab-like spectrum. On the 
other hand, the fluxes of our light curve are in the level of 
~ 100 mCrab at the late phase. This flux level is reasonable 
to detect in ~ 100 s exposure by BAT. Thus, we believe that 
there is a systematic error by u sing the CPL model in the BAT 
flux calculation presented by Campana et al. (2006). We also 
believe that this investigation strengthens our conclusion that 
the Band function is the best represented spectral model of the 
non-thermal emission of GRB 060218. 

2.2. Implications for the nature of the prompt non-thermal 
emission 

2.2.1. Overall spectral properties 



We find that the prompt non-thermal emission of this event 
has spectral properties similar to those of the prompt 7-ray 
emissions of GRBs detected so far. The low- and high-energy 
photon indices as and (3b do not deviate significantly from 
~ —1.0 and ~ -2.5, respectively (Figure[T|). Thes e values are 
quite typical for the prompt emissions of GRBs dPreece et alJ 
2000). Thus the prompt non-thermal emission may be cat- 
egorized into X-ray Flashes (XRFs), which are the transient 
events with properties similar to GRBs exc ept lower spec- 
tral p eak energies and smaller fluences (e.g., Sakam oto et alJ 
2005). The averaged spectral peak energy ~ 4.9 keV and 
the isotropic-equivalent energy # 7l i S0 — 6 x 10 49 erg of the 
non-thermal emission obey the well-known Amati correla- 
tion, which is satisfied by the GRBs and XRFs with known 
redshifts except for outliers such as GRB 980425 and GRB 
031203 dAmati et al.l I2002L 120061) . Moreover, it has been 
reported that this event obeys the lag-l uminosity correlation 
dLiang et al.ll2006dtlGehrels'^aIll2006h . 

We may confirm that the speed of the emitting region was 
close to the light speed. For the emission radius ro, the 
condition of the transparency for the Compton scattering off 
the electrons associated with baryons is given in the non- 
relativistic formulation by 

T — 2 — 7<l, (1) 

2irm p c 2 [3o r 2 

where E KO is the isotropic-equivalent kinetic energy of the 
outflow, oj is the Thomson cross section, and m p is the pro- 
ton mass. Here we assume that the outflow is dominated by 
the kinetic energy, and Eis is larger than the photon energy 
£U i so . Estimating that ro ~ c(3o8t and adopting the values 
£ 7iiso ~ 6 x 10 49 erg and St ~ 3000 s, we obtain f3 a > 0.9. 
This corresponds to the Lorentz factor To > 2. For the 
relativistic wind, the condition of the transparency is given 
by t ^ g7-L isn /(87rmn C 3 r n 3 rn) < 1 dDaigne & Mochkovitchl 
l2002t [Meszaros 2006), where Lj so is the isotropic-equivalent 
kinetic luminosity of the wind. Estimating that L; so > 
Ej.iso/ St and ro ~ To 2 cSt , we obtain no limit on the Lorentz 
factor. (Since the maximum photon energy detected is lower 
than the electron's rest energy m e c 2 , there is no limit on 
the Lorentz factor from the photon annihilation or from the 
Compton scattering off pair-produced electr ons and positrons 
dKaneko et alJl200a lUthwick & Sarill200l .1 Therefore the 
expansion speed of the outflow should be close to the light 
speed, although it is not necessarily ultra-relativistic. 

All these facts show that the prompt non-thermal emission 
of this event is a typical GRB (or XRF), except that its lu- 
minosity is extremely low and it is accompanied by the soft 
X-ray thermal component. This is one of our motivations for 
considering that GRB 060218 arises from a jet like typical 
GRBs in § [3] and [4] In the rest part of this section, we fo- 
cus on the decay phase of the prompt non-thermal emission 
to suggest that its emission process does not stop abruptly but 
decays slowly to connect to the detected anomalous X-ray af- 
terglow. 

2.2.2. The decay phase 

The decay phase of the prompt non-thermal emission is 
well described by F(t) oc f~ 20 and E p (t) oc f" L6 , which im- 
plies a hardness-intensity correlation F oc E p 3 . Such a corre- 
lation, F oc E p ^, can be seen in the decay phases of the GRB 
prompt emissions observed by BATSE, which has a broad 
distribution of the power-law indices £, with 0.6 < C ^ 3 
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dBorgonovo & Rvdd 120011: iRvde & Petrosiar] 120021) . Swift 
satellite has succeeded in observing the decay phases of many 
GRB prompt emissions more deeply than the observations in 
the pre-Swift era. The result is that the distribution of the tem- 
poral ind ices of the flux is al so broad; — 5 < a < — 1, where 
F <xt a (lO'Brien etal.1 120061) . Thus the flux decay of the 
prompt non-thermal emission of GRB 060218 is relatively 
shallow. 

First, we argue that the decay phase is not attributed 
to the kinematical effect due to the curvature of the emit- 
ted region which ceased the emission process suddenly. 
Since the emitting shell should have a curvature, the ob- 
server receives the radiation far from the line of sight af- 
ter the cessation of the emission process. The region at 
higher latitude from the line of sight has a lower veloc- 
ity towards the observer, so that the emission becomes 
dimmer and softer progressively because of the relativistic 
beaming effect. This effect is the so-called curvature ef- 
fect, and have been wid ely studied for the decay phases of 
GRB prompt emissions ( Fenimore et alj 119961: IS ari & Piran 
19971: iKumar & Panaitescul 120001; IRvde & Petrosianl 120021: 
Dermerir2004l) . Swift observations have suggested that the 
decay phases of GRB prompt emissions with the tem- 
poral indices a < -2 are due to the curvature effect 
dLiang et al.ll2006at IZhang et al.ll2006at lLazzati & Begelmanl 
2006l: lYamazaki et al.ll2006r ikobavashi & ZhangJ|2006b. 

Suppose that the shell moving with a bulk Lorentz factor 
Tq ceases the emission at a radius ro and at a time fo and that 
the line of sight is within the shell, i.e., the observer views 
the shell on-axis. The comoving frequency v' is boosted to 
v = v'T> in the observer's frame, where T> = [To(l -/3ocos(9)] _1 
is the Doppler factor. Because the observed time is given by 
t = fo - ro cos #/c, the Doppler factor is related to the observed 
time by 



where t r = to - ro / (c/?o) is the departure time of the shell from 
the central engine. Thus the time-resolved spectral peak en- 
ergy evolves as 

E p <xV<x{t-t r )- 1 . (3) 

The received flux in a given time interval dt is F v dt oc 
j' v ,T> 2 dQ,, where j' v , is the comoving surface brightness, and 
the solid angle of the emitting region is r elated to the observed 
time i nterval by d£l = 2nd(cos9) oc dt (iRvbicki & Lightman 
119791: iGranot et al.lll999t floka & Nakamurall2001l) . Thus the 
received flux is estimated by 

F u oc j' u ,V 2 oc (is') l+ P°V 2 oc v^V^o oc ^(f-f,.)" 1 ^, 

(4) 

where we assume the high-energy Band spectrum f v , oc v' 1+,3b 
because the 15- 150 keV band is above the peak energy E p . 
With the observed value (3b — -2.5, we have F (f) oc (f-f r )" 3,5 
and E p (t) oc (t — t r ) . To reproduce the temporal index ob- 
tained from our analysis, t r < is required for F(t), while 
in contrast, t r > is required for E p (t). Therefore F (f) and 
E p (t) cannot be fitted simultaneously by this model. One can 
also consider the structured jet model in w hich the spectrum 
and brightness vary in the angular directi on (Rossi et al.l2002t 
Zhang & Meszarosl 120021: lYamazaki etail 120061: iDvks et all 
2006). However, the shallower flux decay requires a jet with 
brighter rim, which seems implausible. 

Thus the emission process is attributed not to the curva- 
ture effect of the shell which ceases the emission abruptly, 



but to slowly decaying emission processes. The first possi- 
bility of the decaying e mission is the synchrotron emission 
from the external shock (Rees & Meszaros 1992; Sari et al. 
1998). Since the high-energy range of the spectrum will be 
above the cooling frequency, the high-energy photon index 
Pb — -2.5 corresponds to the energy distribution index of the 
accelerated electrons p ~ 3. Then the flux should decay as 
F(t) oc t~-~'* p ~ f" 1,75 . The characteristic frequency v,„ evolves 
as f" 3 / 2 . These temporal behaviors are not inconsistent with 
the results of our analysis. However, the standard model gives 
much lower characteristic frequency v m than the observed one 
~ 10 18 Hz at the deceleration time t ~ 500 s with the X-ray 
luminosity L x ~ 10 47 erg s" 1 . Therefore this model is not so 
appealing. 

The second possibility is that the decaying emission is at- 
tributed to the central engine activity. The decaying flux in 
15- 150 keV range, if extrapolated as F(t) oc ?~ 20 , becomes 
~ 10~ n erg s" 1 at t ~ 10 4 s. This flux is comparable to that of 
the detected anomalous X-ray afterglow in 0.3- lOkeV range, 
decaying as Fx oc f~ with the photon index 0x — -3.2, 
which cannot b e explained within the external shock syn - 
chrotron model (Campana et al. 2006; Soderberg et al. 2006). 
If the anomalous X-ray afterglow is a continuation of the 
prompt non-thermal emission (i.e., the X-ray emission de- 
tected at t > 10 4 s is not an afterglow but the decaying prompt 
emission), the prompt non-thermal emission of this burst is 
very long (> 10 6 s) and it is necessary that the photon in- 
dex var ies from /3g — -2.5 to fix — -3.2. G hisellini et alj 
(2006b) have also suggested such a scenario and shown that 
a synchrotron inverse-Compton model could reproduce the 
UV/optical thermal component as well as the prompt non- 
thermal emission and the anomalous X-ray afterglow. Here 
we note that the flux decays steeper than t~ l after the peak 
time t ~ 500 s, so that the time-integrated radiation energy 
does not increase so much after the peak time. Thus the total 
radiation energy of the prompt non-thermal emission remains 
£ 7 , iso ~ 6 x 10 49 erg. 

Therefore the engine of this event could be active for at 
least 10 6 s. Recent detailed observations of GRB X-ray af- 
terglows have revealed that a number of bursts show erratic 
X-ray flares, which may be linked to the restarting s of the cen- 
tral en gine after the prompt emission phase (e.g . JZhang et alj 
2006l lBurrows et al.ll2005bl: lloka et al.1 12005). Since some 
X-ray flares are discovered ~ 10 5 s after the trigger, the en- 
gines of typical GRBs could be also active for > 10 5 s. How- 
ever, they are spiky and intermittent activities, while this event 
shows the smooth and power-law activity. We speculate that 
this indicates that the engine of this event is different from 
those of typical GRBs which are believed to be black holes. 
Although this is just a speculation, it agrees with the analy- 
sis of the SN component of this event, which suggests that 
the progenitor star was less ma ssive and collapsed to a neu- 
tron star instead of a black hole (Mazzali et al. 2006). Further 
discussion will be presented in §[5j 

3. A JET MODEL OF THE RADIO AFTERGLOW 

The afterglow of this burst shows chromatic light curves. 
The UV/optical afterglow is quite dim and d ominated by the 
thermal component and the S N component (ICampana et alJ 
120061: iGhisellini etail l2006bl) . The X-ray afterglow has a 
spectrum much steeper than tho s e of typical GRB X-ray af- 
terglows (ISoderb erg et al. 2006; Butler 120061) . and it could 
be a continuation of the prompt non-thermal emission (see 
§ 12.2.2b . Only the radio afterglow seems rather typical and 
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can be explained by the standard external shoc k synchrotron 
model (ISoderberg et al.ll2006t iFan et al J 120061) . In this sec- 
tion, we show that the available radio data may arise from an 
initially collimated outflow. In the spherical outflow model, 
it seems difficult for a SN shock to provide the relativistic 
spherical ejecta with sufficient energy, as stated in § [TJ In the 
jet model, taking into account the UV/optical and X-ray after- 
glows, we will constrain the initial opening angle and Lorentz 
factor of the jet as do — 0.3 and r ~ 5, respectively. We also 
show that the off-axis scenario in the jet model is less likely 
because it requires higher 7-ray efficiency than the on-axis 
scenario. 

3.1. The radio afterglow 

For 2-22 days the afterglow flux at 8 .46 GHz is well de- 
tected and shows the power-law decay as t~~° 85 . The spectrum 
at t ~ 5 days has a peak between 1 .43 GHz and 4.86 GHz, and 
this peak can be interpreted as the self-absorption frequency 
v a in the standard model. This peak would not be the typi- 
cal synchrotron frequency v m because of the small 1 .43 GHz 
flux. Thus the 8.46 GHz is most likely in the frequency range 
v m < v a < v < v c , where v c is the cooling frequency. 

In the standard model, the temporal index of the flux at 
each frequency range can be calculated for several cases; 
for a relativistic/non-relativistic blastwave expanding into cir- 
cumburst medium of constant/wind-profile density, or for a 
sideways expanding relativistic blastwave. The temporal in- 
dex of the flux at the frequency range v m < v a < v < v c for 
each case is — §/> + § (relativistic blastwave, constant density), 
— f/?+fg (non-relativistic blastwave, constant density), —\p + 
i (relativistic blastwave, wind-profile density), -gp+f (non- 
relativistic blastwave, wind-profile density), and — p (sideways 
expanding), where p is the index of the e nergy distribution 
funct i on of the accelerated electrons (e.g., [L ivio & Waxman 
2000 [lSari et al.ll 19981 fl999t IChevalier & Lill 1 999t iFraifetaT 
2000). For each case to reproduce the observed temporal in- 
dex -0.85, p is required to be ~ 2.1,2.0,1.5,1.4, and 0.9, 
respectively . The values p > 2 are typical for the GRB after- 
glows fe.g. JPanaitescu & Kuma^2002t lYost et al.ll2003l) . and 
for this reason, we focus on the first two possibilities (i.e., the 
constant density cases). Note, however, that the possibilities 
p < 2 (i.e., the wind-profile density cases and the sideways 
expanding case) cannot be ruled out. 

For the two possibilities of the standard model, the available 
radio data requires the following three conditions: 

(a) The spectrum peaks at v a ~ 4 x 10 9 Hz at t ~ 5 days. 

(b) The flux at 22.5 GHz is ~ 0.25 mJy at t ~ 3 days. 

(c) The cooling frequency is v c < 5 x 10 15 Hz so that the 
synchrotron spectrum of the external shock electrons 
does not dominate the detected anomalous X-ray after- 
glow. 

3.1.1. Relativistic blastwave model 

IFan et al.l d2006l) have modeled the first possibility, i.e., 
the synchrotron spectrum from the relativistic blastwave with 
p ~ 2.1 expanding into the constant density medium, and de- 
rived three constraints on the physical parameters of the af- 
terglow from the above three conditions (a) - (c) (their Equa- 
tions (8) -(10)). The satisfying parameters they have sug- 
gested are the isotropic-equivalent kinetic energy of the out- 
flow Ekjao ~ 10 50 erg, the number density of the circumburst 



medium n ~ 10 2 cm" 3 , and the ratios of the accelerated elec- 
trons energy and the magnetic energy to the shocked thermal 
energy e e ~ 10~ 2 and ~ 10~ 3 , respectively. With these pa- 
rameters, the Lorentz factor of the outflow is estimated as 
r»2(f/l day)" 3 / 8 . If the opening angle of the outflow is 
9q< 1 , the outflow will begin to exp and sideways a t t < 6 days 
and the radio light curves will break (Rhoads 1999). However, 
the actual light curve does not show such a break, so that it is 
concluded that 6q > 1 in this model, i.e., the outflow is spher- 
ical. In this case, the outflow becomes non-relativistic at t ~ 6 
days, and the radio flux varies from f -085 to t~ l A . This is co n- 
sistent with the available radio light curves (IFan et alf2 006). 



The parameters E^\ 



10 48 erg, n ~ W cm" 3 , e e 



10 , and eg ~ 10 _1 also satisfy the above thre e constraints 
and fit the radio data dSoderber g et al.l 120061) . However, 
since the isotropic-equivalent 7-ray energy is £ 7 ,iso — 6 x 
10 49 erg, the inferred 7-ray efficiency is extremely high; 
?7 7 = -E-y.iso/ (Eyjso+Ekjuo) ~ 98%, and thus this parameter set 
is implausible. 

The relativistic spherical outflow is not consistent with 
the standard collapsar scenario, in which a relativistic col- 
limated outflow digs a narrow hole in the progenitor star 
without loading so much s tellar matter dWooslevI 119931 ; 
iMacFadyen & WoosleyH l999). Thus one may consider a sce- 
nario in which a spherical outflow is accelerated to a mildly 
relativistic speed as the SN shoc k sweeps the stellar enve- 
lope with steep density gradient dMatzner & McKed 119991 : 
iTan et al]l2001h . However, it seems impossible that the re- 
quired kinetic energy E% t i so ~ 10 50 erg is transferred into the 
relativistic ejecta by the SN w ith the total kinetic energy 
~ 10 51 erg dMazzali et al.ll200"6l) . as discussed in §Q] For this 
reason, we conclude that the first possibility is quite unlikely. 



3.1.2. Non-relativistic blastwave model 

Now we investigate the second possibility; the radio af- 
terglow may arise from the non-relativistic phase of an ini- 
tially jetted outflow with p ~ 2.0 expanding into the con- 
stant density medium. The relativistic jet is decelerated 
and subsequently expands sideways, and finally becomes 
non-relativistic at a certain transition time t s . The transi- 
tion time is estimated by requiring that the speed of the 
Sedov-Taylor blastwave is close to the light speed, i.e., (3 = 
■^[Eit/inmnt 3 )] 1 / 5 ~ 1, where Ei is the kinetic energy of the 
spherical blastwave (see Livio & Waxman 2000): 



:7.5xl0 6 s'^' 51 



1/3 



(5) 



Here (and hereafter) we have adopted the notation Q = 10 X Q X 
in cgs units. Since the available 8.46 GHz lightcurve does not 
show the sideways expanding relativistic phase t~ p ~ f" 2 , we 
require t , < 2 days. 

Three constraints on the physical parameters of the after- 
glow are given by the conditions (a) - (c) in the first part of 
this subsection (§ 13. U . Assuming that the minimum Lorentz 
factor of the accelerated electrons is j m ~ 10~ l f3 2 e e m p /m e , 
we may derive the following three constraints on the physical 
parameters: 



'/ 3 £ l/V/3 



1. 



0.75 F 1.3-0.45 o v i n - 
-16R-2 Eb51 n ~ ' x 1U 



n 

0.6„0.9 



>0.4. 



(6) 
(7) 
(8) 
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The results are similar to th ose derived in t he relativistic 
blastwave model (see § 13.1.11 iFan et al.ll2006l) . These con- 
straints and the requirement of t s < 2 x 10 5 s are satisfied by 
£ t ~2x 10 48 erg,« ~ 10 2 cm" 3 ,£ e ~ KT 1 , and e B ~ 10" 1 . 

Then the isotropic-equivalent kinetic energy before side- 
ways expansion is estimated by E^ = 2E k j0Q ~ 4 x 

10 48 9o 2 erg. To obtain a rea sonable 7-ray efficiency 
W-y = E~y.j s0 /(E^ so + E ti*r.) < 0.5 dLlov d-Ronning & Zhang 
2004tlloka et al. 2006; Fa n & Piranlf2006l: iGranot et al.ll2006t 
Zhang^eLaLl|2006b), the opening angle 0q J$ 0.3 is favorable. 
In the following section, we will further constrain the opening 
angle and also the Lorentz factor of the jet using the observed 
UV/optical and X-ray emission. 



jet is determined by 



1/2 



L. 



(9) 



When t s ~ 2 x 10 5 s and tj ~ 2 x 10 4 s are adopted, we ob- 
tain #0 ~ 0.3. If 8o is smaller, the ratio t s /tj is larger. At the 
Blandford-McKee phase, the Lorentz factor of the jet evolves 
as r oc f" 3 / 8 , and finally becomes Y ~ f?o~ at the jet-break 
time tj. Thus the initial Lorentz factor of the jet is determined 
by 



^dec 



-3/8 



(10) 



3.2. Constraints from the UV/optical and X-ray afterglow 

In our model (£ t ~ 2 x 10 48 erg,n ~ 10 2 cm" 3 ,£ e ~ 10" 1 , 
and eg ~ 10 _1 ), the characteristic parameters at t ~ 5 days 
are calculated as iv max ~ 2 x 10~ 25 erg cm" 2 s -1 Hz -1 , 
i/„~5x 10 6 Hz, v a ~ 4 x 10 9 Hz, and z/ e ~ 7 x 
10 13 Hz. Thus the optical ^/v flux at this time is ~ 
l0 l5 F u . m ^c/is m y - 5 (l0 l5 /v c y l ~ 4 x 10" 15 erg cm" 2 s" 1 . If 
we trace back to earlier times through the standard jet model, 
the optical flux gets much higher than this value (see Fig- 
ure The detected optical flux at 10 4 - 10 5 s is dominated 
by the thermal component with vF v ~ 10~ n erg cm" 2 s" 1 
(Ghisellini et al. 20063 ICampana et alll2006l) . The require- 
ment that the synchrotron flux from our jet should not exceed 
this flux gives constraints on the initial opening angle 0q and 
the initial Lorentz factor To of the jet. 

The external shock of the jet progressively experiences the 
Blandford-McKee evolution, the sideways expanding evo- 
lution, and the Sedov-Taylor evolution. For each evolu- 
tion phase the characteristic frequencies of the synchrotron 
spectrum with p ~ 2.0 vary as v m oc f~ L5 ,f~ 2 , and f" 3 ; 
v a (> v m ) oc r 2 / 3 ,r ' 5 , and r 2 / 3 ; v„ (< v m ) oc f°,r a2 , and 
t L2 ; and v c oc r°- 5 ,f° and r 02 (e.g.. ISari et aDfT998l fl99l 
iLivio & Wa xman 2000). If we trace back to t ~ 10 4 s, v m ,u a , 
and v c all remain lower than the optical band. Thus for 
t > 10 4 s, the optical band is at the range v > v c , and the opti- 
cal energy flux evolves as t~ l at the Blandford-McKee phase, 
f" 2 at the sideways expanding phase, and t~ l at the Sedov- 
Taylor phase. With these behaviors of the optical flux from 
the external shock, we can suggest a following possible jet 
scenario (see Figure [3] for a schematic picture of the optical 
lightcurve). The external shock of the jet is decelerated at 
fdec — 6 x 10 3 s, begins to expand sideways at tj ~ 2 x 10 4 s, 
and shifts to the Sedov-Taylor expansion at t s ~2x 10 5 s. The 
peak optical flux is vF v ~ 3 x 10~ 12 , which is less than the ob- 
served one, so that the optical afterglow estimated within our 
jet model is not inconsistent with the available optical data. 

Since the X-ray band is also above v c and p ~ 2.0, the 
X-ray vF v flux is comparable to the optical one. Thus this 
external shock emission does not over whelm the observed 
ano malous X-ray af t erglo w shown in Cam oana et al.1 (2006) 
and lSoderberg et al.l (120061) . 

Now we give the reasonable values of 0q and To. At 
the sideways expansion phase, the opening angle of the jet 
evolves as = Y' 1 oc t l l 2 (lRhoadslll999HSari et al.lfT999l) . and 
finally becomes ~ 1 at the transition time t„ to the Sedov- 
Taylor expansion phase. Thus the initial opening angle of the 



When tj ~ 2 X 10 s and fdec — 6 x 10 s are adopted, we ob- 
tain To ~ 5. If To is larger, the ratio f//fdec is larger. The 
upper bound of the optical data does not allow the values of 
ts/tj and tj/tdec larger than those we have adopted, and there- 
fore it does not allow the smaller 0q and the larger To. Since 
the consideration about the 7-ray efficiency requires 0q < 0.3 
(see § 13. 1.2b . the opening angle is restricted to a narrow range 
around 6*o — 0.3. 

In a summary, we have obtained a jet model for the radio 
afterglow of GRB 060218, with O ~ 0.3, T ~ 5, E kMo ~ 
4 x 10 49 erg, and r; 7 ~ 0.6, which is compatible with the 
UV/optical and X-ray data. The collimation-corrected en- 
ergy of the jet is Ej ~ (^ 4so + £ , ^, iso )6'o 2 /4 ~ 10 48 erg. As- 
suming that the duration of the prompt non-thermal emission 
St ~ 10 3 s is the active time of the central engine ST, the lu- 
minosity of the jet is estimated by Lj ~ EjjST ~ 10 45 erg s -1 . 

3.3. Off-axis scenario 

In the above arguments in §[2]and[3j we have assumed that 
the observer views the jet from the on-axis direction and the 
luminosity of the prompt non-thermal emission is intrinsically 
low. It is possible that the low luminosity of GRB 060218 is 
attributed to a typical bright GRB jet viewed off-axis (e.g., 
Yama zaki et alj|2003t iRamirez-Ruiz 1 et alj|2005t IGranot et all 
2005b iToma et alJl2005h . iowever, we find that the off-axis 
scenario is less likely than the on-axis scenario since higher 
7-ray efficiency is required. 

Since the outflow is spherical when the radio afterglow is 
observed, the off-axis scenario is the same as the on-axis one 
in this phase. Thus the isotropic-equivalent kinetic energy is 
determined by E k>iso = 2E k /0 Q 2 ~ 4 x 10 48 Q ~ 2 erg. In the off- 
axis scenario, the isotropic-equivalent 7-ray energy Eyj so on 
measured if viewed on-axis is larger than the observed one 
E 7i i S0 ~6x 10 49 erg. So the 7-ray efficiency is expected to be 
higher than that estimated within the on-axis scenario (e.g., 
IToma et al.f2006l) . 

To obtain the smaller 7-ray efficiency, we could consider 
0{) < 0.3. However, the smaller 0o leads to even higher 7- 
ray efficiency, as we will see below. The smaller #0 requires 
tj < 10 4 s, and then the optical flux, if received from the on- 
axis direction, exceeds the constraint from the detected optical 
afterglow (see Figure [3). To escape the constraint, the syn- 
chrotron flux from the off-axis jet should peak at some time 
t v after tj. The peak occurs when the line of sight enters the 
relativistic beaming cone of the emission from the edge of the 
sideways expandingjet, i.e., V ~ 6*+r _1 ~ 2T" 1 . The Lorentz 
factor varies as T oc t~ l l 2 in the sideways expansion phase and 
becomes Y ~ 1 at t s , and thus the following equation is satis- 



A collapsar jet model of GRB 060218 



7 



fied: 



1/2 



1. 



(ID 



Since t v should be greater than 10 4 s, we obtain 9 V > 0.4. 
The isotropic-equivalent 7-ray energy measured if viewed on- 



axis is roughly estimated by 



[To(9 v -9o)] E~ 



> 

10 47 r 6 erg, where we have used 9 V > O . With r > 9 ~\ 
we obtain £ 7i i S o,on > 10 47 #o~ 6 erg. For 6q < 0.3, this leads 
to higher 7-ray efficiency than the on-axis scenario; rj 7 > 0.7. 
Therefore we may conclude that the off-axis scenario is un- 
likely for this event. 

4. A COLLAPSAR MODEL OF GRB 060218 

In the previous section, we have obtained a possible jet 
model that is consistent with the available afterglow data. Our 
jet model requires an opening angle 9q ~ 0.3, which is some- 
wh at larger than those of typical cosmological GRBs (e.g., 
see iGhirlanda et al.l l2004). and an extremely small luminos- 
ity Lj ~ 10 45 erg s _1 , which is about 10 5 times smaller than 
those of typical cosmological GRBs. Then one may won- 
der whether such a wide and low-luminosity jet can pene- 
trate a progenitor star. It is possible that a jet becomes non- 
relativis tic if it loads so much stellar matter with a large cross 
section. iMatzerl (120031) has discussed analytically several con- 
ditions for making a hole in the collapsar model. We extend 
his theoretical considerations to conclude that an adiabatic 
cold jet is excluded, but a non-adiabatic hot jet is appropriate 
for this event. (Our argument can be also applied to typical 
GRB jets, and some wide GRB jets also would need to be 
hot.) 

4.1. The motion of a jet head 

The standard collapsar model assumes that a black hole or a 
neutron star with an accretion disk is formed after an iron core 
of the massive progenitor star collapses and that the system 
produces a jet dWooslevll993TlMacFadven & Woosievll 19991) . 
Since the free-fall timescale of a stellar envelope is longer 
than, or comparable with, the crossing timescale for the jet to 
propagate within the progenitor star and to hit its surface, the 
stellar envelope would disturb the progress of the jet. 

Consider the progress of a relativistic jet outward through 
the stellar envelope. Two shocks form: a reverse shock re- 
ducing the jet speed and increasing its internal energy; and 
a forward shock that propagates into the surrounding stellar 
envelope giving it internal energy. As a result, there are four 
distinct regions in this system: the propagating jet; the head 
of the jet lying between the two shocks; the stellar envelope; 
and a cocoon consisting of shocked je t and shocked ambient 
material (see Figure 2 of lMatzerl2003l). This system is simi lar 
to classical double radio sources (Bege lman & Ciofrj l989). 

Two conditions are required for a jet to break out the star 
relativistically. First, the Lorentz factor of the jet head T/, 
should be smaller than the inverse of the opening angle of 
the jet 9. If so, the shocked material may escape sideways 
to form the cocoon, and the jet may avoid baryon loading. 
Second, the speed of the jet head should be larger than that of 
the expanding cocoon, 

Ph > (12) 

If this condition is violated, the cocoon expands spherically 
around the jet and finally explode the star, producing a non- 
relativistic dense spherical outflow. When the above two con- 
ditions are satisfied, we may consider the longitudinal balance 



of the momentum flux in the frame of the jet head, 

w;r/r A 2 (& - I3 h f + Pj = W a T h 2 [3 h 2 + Pa . 



(13) 



Here the subscripts j and h describe the jet and the jet head, 
respectively, and w(= e + p), p, and e are the enthalpy, the 
pressure, and the energy density including rest energy density 
pc 2 , respectively. For stellar envelopes, p a <C p a c 2 and w a = 
p a c 2 are good approximations. 
If r ; is significantly larger than unity, pj can be neglected 

in the left-hand side of Equation dTJl i. and we obtain (Pi,~ l - 
l)" 2 = wjT 2 /(p a c 2 ) = Li so /(4nr 2 p a c 3 ), where at the last equal- 
ity we introduce the isotropic-equivalent luminosity L; so of the 
jet. For the jet to drive the explosion of the star, this equation 
should be satisfied when r = R, where R is the radius of the 
progenitor star, so that (P h ~ l - 1)" 2 w l0~ 3 L iso , 5 i(R/R Q )~ 2 p a ~ l . 
This factor is much less than unity for the parameters of typ- 
ical GRBs (and so for the low-luminosity GRBs), so that 
the jet head is found to be non-relativistic. Thus, with the 
collimation-corrected luminosity of the jet Lj = WjT 2 cirr 2 9 2 , 
we obtain 



Ph 



Trr 2 9 2 p a c- 



L/2 



(14) 



We take p a ~ M /(4irR 3 /3) for simplicity, where M is the total 
mass of the progenitor star. 

4.2. Cocoon structure 

The cocoon drives a shock into the stellar envelope, which 
expands non-relativistically at the velocity j3 c given by the 
transverse balance of the cocoon pressure and the ram pres- 
sure of the stellar envelope, 



p c = p a C P c 



(15) 



The pressure of the cocoon p c is roughly constant away 
from the jet head, and is approximated by £ , ; n /(3V c ) with 
the thermal energy E m deposited in the cocoon and the co- 
coon volume V c . Since the jet head is non-relativistic, most 
of the kinetic energy of the jet gets thermalized through 
the reverse shock. The cocoon energy is the energy caught 
by the jet head, which is calculated by E m = JqLj-^- ~ 
Ljjp-. The cocoon length is given by r ~ cf3ht as long 
as fih > Pc is satisfied, and the cocoon width is given by 
R c ~ cft c t. So the cocoon volume is roughly V c — (n/3)R c 2 r ~ 
(Tr/3)r 3 p 2 /p h 2 . Therefore the cocoon pressure can be written 
by p c ~ LjPh/(irr 2 cPc 2 )- Equation < fT3T > gives us the following 
equations: 



Pc'- 



irr 2 p a c 2 



Ph 



-^PaCPh 



1/4 



1/2 



(16) 



(17) 



which describe the velocity and the pressure of the cocoon 
with the parameters of the jet and the stellar envelope within 
the assumption of Ph > P c - 

4.3. Cold jet 

Here we examine an adiabatic jet that propagates ballisti- 
cally in the progenitor star (see also Meszar os & Reesll2001t 
IWaxman & Meszarosll2003l) . An important point is that the 
opening angle and the Lorentz factor of the cold ballistic jet 
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do not change after exiting the progenitor star. If the open- 
ing angle and Lorentz factor of the jet at the breakout are 6*b r 
and Tbr, they are also 9q = #br and Tq = after exiting the 
star. In contrast, a hot jet which propagates non-adiabatically 
by the interaction between the jet and the cocoon and keeps 
dominated by the internal energy changes its opening angle 
and Lorentz factor through the free expansion outside the star. 
The hot jet will be discussed in the next subsection. 

For the jet to drive an explosion, the consistency f3f, > f3 c 
should be satisfied. From Equations (fl4] i and (fTol l. this leads 
to the constraint on the opening angle of the jet, 



< 



irr 2 p a c 3 



1/6 



(18) 



For GRB 060218, the underlying SN is Type Ic, which 
implie s that the progeni tor star is a C/O Wolf-Rayet star 
(e.g., iMazzali et al.ll2.Q06b . Its radius and mass are typically 
R ~ 10 11 cm and M ~ 10 34 g, respectively. Then the av- 
eraged mass density p a is ~ 1 g cm" 3 . For the luminosity 
L;~ 10 45 ergs -1 , the opening angle of the jet at the breakout 
(i.e., when r = R) is constrained to 6^ < 0.03, so that 9q < 0.03 
within the cold jet scenario. The value of the opening angle 
suggested in our jet model (Oq ~ 0.3) violates this constraint. 
Therefore the cold jet scenario is excluded for this event. 

A wide cold jet may be allowed if the jet is being launched 
long after the cocoon explodes the star, since the ambient den- 
sity p a drops through the expansion in Equation d 1 8t . How- 
ever the star must expand from R c ~ 10 11 cm to ~ 10 17 cm in 
order that 9q ~ 0.3 is allowed, which is unlikely. 

4.4. Hot jet 

Next we consider the possibility that the jet is dominated 
by the internal energ y throughout the propagatio n in the pro- 
genitor star (see also lLazzati & Be gelman 2005). At the jet- 
cocoon boundary, oblique shocks and shear instabilities may 
occur and dissipate the kinetic energy of the jet. This sit- 
uation can be se en in several numerical simulations of the 
collapsar model d Alovetai] |2000t IZhang et all 120031 l2004t 
lUmeda et al.ll2005UMizuta et al.ll2006l) . 

In contrast to the cold jet, the hot jet expands freely after 
exiting the star. As a result, the opening angle is determined 
by the Lorentz factor r ; ex just before exiting the star. In 
the comoving frame of the hot material moving outward rela- 
tivistically with the Lorentz factor r ; ex , the hot material will 
expand freely and finally get the Lorentz factor V = e/pjC 2 . 
In the laboratory frame, the material will be beamed into the 
opening angle 

do ~r Jja r l , (i9) 

and the Lorentz factor of the material is given by To ~ 
r'r ;> x(l + /3'/3>.ex) ~ 2rT;, ex . Note that T O O ~ 2T' > 1 is 
satisfied. 

For GRB 060218, Tj ex ~ 6~q ~ 3 is required to reproduce 
9o — 0.3. For the initial Lorentz factor of the jet exiting the 
star to be To ~ 5, V = e/pjc 2 ~ 1 - 2 is favorable. Therefore 
the jet of this burst may originate from a mildly hot jet in 
which the internal energy e is comparable with the rest energy 
PjC 2 before exiting the star. 

4.4.1. The breakout timescale 

Let us estimate the breakout timescale. We expect that the 
breakout timescale is not much larger than the active time of 
the central engine after the breakout ST ~ 10 3 s. Otherwise, 



the engine has to stop suddenly just after breaking out the star, 
and this seems unlikely. 

For the hot jet, pj ^> pjc 2 and Wj = 4pj are good approxi- 
mations, and thus pj = Lj/ (4T j 2 or r 2 9 2 ). Before the breakout, 
the opening angle of the hot jet is determined by the transverse 
pressure balance between the jet and the cocoon: 



Pi= Pa- 
using Equations ( TT4T > and (fTTT i. we obtain 



256irr 2 p a c- 



1/6 



r - 4/3 



(20) 



(21) 



which satisfies the consistency of (3f, > (3 C (Equation (fT8l). 
For GRB 060218, adopting the values Lj ~ 10 45 erg s" 1 , R ~ 
10 11 cm, p a ~ 1 g cm" 3 , and Tj ~ 3, we obtain the opening 
angle at the breakout time as 6*b r ~ 3 x 10~ 3 . 

The velocities of the jet head (Equation ( TBI ) and the co- 
coon (Equation ( fT6l l) are calculated by using Equation < f2~TT >. 



\6Lj V /3 t^.4/3 




(22) 
(23) 



The ratio of the two velocities is /3/,//? c = 2Tj. Then, if Tj is 
close to unity, we have (3 C ~ Pi, and hence the cocoon ex- 
pands quasi-spherically. This situation is in a good agree- 
m ent with recent seri es of numerical simulations performed 
by iMizuta et al.l (120061) . They have argued that the morphol- 
ogy of the explosion depends on the Lorentz factor Tj o of the 
jet given at the inner boundary: when the Lorentz factor is 
high (r^-.o ^ 3), the high pressure cocoon collimates the out- 
flow to form a narrow, relativistic jet; and when the Lorentz 
factor is low, on the contrary, the outflow is not collimated and 
expands quasi-spherically. 

The breakout timescale is determined by tbi = R/ (cf3h.br), 
where j3h,bt is the velocity of the jet head (Equation d22l) ) ob- 
tained at r = R: 



tbr 



-1/3 



^5/3 r .-4/3^ 



(24) 



For GRB 060218, we obtain fb r ~ 300 s, which is a reasonable 
value. The hot jet of this event takes fbr ~ 300 s to break out 
the progenitor star and lasts further ST ~ 10 3 s to eject the 
material producing the prompt emission. 

4.4.2. A possible additional widening effect 

After the breakout time, the opening angle of the jet be- 
fore exiting the star might become wider than the value 0b T ~ 
3 x 10~ 3 , which is determined by the pressure balance be- 
tween the jet and cocoon. The cocoon expands freely outward 
from the star after the breakout time, and finally the cocoon 
pressure becomes negligible in R/(c/^/3) ~ 10 s. The cocoon 
shock that will sweep the stellar envelope is not so strong, and 
the envelope would remain cold. After the cocoon disappears, 
the jet opening angle would be determined by the transverse 
balance between the jet pressure and the ram pressure of the 
stellar envelope matter. The balance equation in the comoving 
frame of the propagating jet is 



p j = (T j p a )(r9f 



(25) 
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where r is the radial coordinate in the laboratory frame. A 
Lorentz factor appearing in the right-hand side is due to the 
Lorentz transformation. Substituting pj = Lj/{4Y 2 irr 2 9 2 c), 
we find that 9 gradually increases since the jet pushes the stel- 
lar envelope. Let Lj and Tj be constant and integrate over 
time, and we obtain 

/ T- \ 1/4 

0=1 — ) r- l Yf^8T 1 ' 2 . (26) 
\irp a cj 

For GRB 060218, if we adopt the values Lj ~ 10 45 erg s _1 , 
p a ~ 1 g cm" 3 , r = R ~ 10 11 cm, r, ~ 3, and ST ~ 10 3 s, we 
obtain 9 ~ 0.1. This is less than the opening angle after free 
expansion 9q — 0.3 in Equation ( fT9] l, so that the widening ef- 
fect is not important in this event. In some cases the widening 
effect could dominate the free expansion effect. 

5. SUMMARY AND DISCUSSION 

We have investigated whether GRB 060218 arises from a 
collimated jet. So far the lack of the jet break has led to the 
interpretation that the outflow of this event is spherical, and 
thereby the outflow is not the standard collapsar jet but the 
outermost parts of the stellar envelope that the SN shock ac- 
celerates to a m ildly relativistic speed ( Soderberg et alJl2006t 
iFan et al.ll2006l) . However, we have shown that the avail- 
able radio data may be interpreted as a non-relativistic phase 
of an initially collimated outflow within the standard exter- 
nal shock synchrotron model, and that the jet model with an 
initial opening angle 9q ~ 0.3, Lorentz factor To ~ 5, and a 
collimation-corrected luminosity Lj ~ 10 45 erg can explain 
the radio data and is compatible with the UV/optical and X- 
ray data. This model is more natural than the initially spheri- 
cal outflow model, because in the latter model, the relativistic 
ejecta for the prompt and afterglow emission with ~ 10 50 erg 
could not be produced b y the underlying SN of the total ki- 
netic energy ~ 10 51 erg (Maz zali et al.ll200"6l) . Furthermore, 
the jet model is supported by the recent report of the detec- 
tion of optical linear polariza tion in the SN component of this 
event (Gorosab el et al.l 120061) . We also show that the jet of 
this event can penetrate the progenitor star by extending the 
analytical considerations by Matzej] d2003l) . The jet would 
be relativistically hot in the progenitor star and hence expand 
freely outside the star into the relatively wide opening angle. 
The off-axis scenario within the jet model is unlikely because 
it requires an unrealistically high 7-ray efficiency. 

With the jet opening angle 9q ~ 0.3, the collimation- 
corrected 7-ray energy of this event is E 1 ~ 3 x 10 48 erg, and 
it makes this event consistent with the Ghirla nda correlation 
dGhirlanda et alJl2004tlGhisellini et alJl2006ah . 

For the prompt emission, we have analyzed the BAT and 
XRT data and found that the non-thermal component of the 
prompt emission may be fitted by the Band function similarly 
to other GRBs. The low- and high-energy photon indices are 
almost constant with as — -1 and (3b — -2.5, respectively, 
which are quite typical features of GRBs. The 15- 150 keV 
flux of the prompt non-thermal emission shows a relatively 
shallow decay, and we show that this decay is not due to the 
curvature effect of the emitting shell that ceases the emission 
process suddenly. The decay of the non-thermal emission may 
directly connect to the anomalous X-ray afterglow detected 
up to t ~ 10 6 s. If this is correct, the central engine might 
be active for > 10 6 s. In summary, the low-luminosity GRB 
060218 has a typical prompt non-thermal emission and may 
originate from a standard collapsar jet possibly driven by a 
long-acting central engine. 



A number of Swift GRBs show late-time X-ray flares, which 
may be attributed to the long intermittent activities of t heir 
central engines (e.g JZhang etalll2006atlfoka et alj|2005l) . In 
contrast, the engine of this burst might has the long power- 
law activity. Then it is possible that the engine of this burst is 
different fro m those of typi c al GR Bs which are believed to be 
black holes. Mazza lFet al.l d2006i) have performed a detailed 
modeling of the spectra and light curve of the SN compo- 
nent and argued that the progenitor star of this event had a 
smaller mass than other GRB-SNe, suggesting that a neutron 
star rather than a black hole was formed as the central engine 
of the jet. The intrinsic rate of such low-luminosity GRBs 
would be larger than the local rate of typical cosmologi- 
cal GRBs {S oderberg et a fl2006t IPian et alj2006t ICobb et al.l 
l2006t iLiang et al.ll2006bl) . For these reasons, low-luminosity 
GRBs might be a distinct GRB population involving neu- 
tron star engines. We speculate that massive progenitor stars 
form black holes at the core collapse which produce highly- 
relativistic jets making high-luminosity GRBs with strong 
spiky prompt emissions and flares, while less massive progen- 
itor stars form neutron stars which produce mildly-relativistic 
jets making low-luminosity GRBs with weak smooth prompt 
emissions. 

If the opening angles of the low-luminosity GRBs are 
around ~ 0.3, the true rate of the low-luminosity GRBs with 
a beaming correction is ~ 10 3 Gpc~ 3 yr" 1 . The local rate 
of Type Ibc SNe is R SN ~ 10 4 Gpc" 3 vr"' dSoderberg et al l 
l2006t Dahlenetllll2004l rCappellaro et al.lll999l) . Then the 
low-luminosity GRBs might be created at ~ 10% rate of Type 
Ibc SNe. By comparison, the collimation-corrected rate of 
the typical cosmological GRBs is ~ 1 % of that of Type Ibc 
SNe. Note that it is also suggested that the birthrate of Gala c- 
tic magnetars is ~ 10% of SN rate dKouveliotou et all 1998). 

Now if the neutron star loses its rotational energy mainly 
through magnetic dipole radiation, and a fraction of energy 
/ is transferred t o the jet, then the luminos ity of the jet is 
approximated by (Shap iro & Teukolskvl 1983b 

L ; ~/Lo(l + ^) 2 , (27) 

where 

L = ^f^L ~ 10 47 erg s" 1 BA/Po,^, (28) 
6c J 

t = ~* C * 2 ~ 10 3 s h 5 Bif 2 R S}6 -*P % - 2 2 . (29) 
B £ R S S io 

Here B is the magnetic field strength at the pole, R s is the 
neutron star radius, Slo is the initial angular frequency, Pq = 
Itt/VLq is the initial rotation period, and / is the moment 
of inertia of the neutron star. The characteristic spin-down 
timescale t might be the peak time t ~ 10 3 s of the prompt 
non-thermal emission, and the temporal decay of the spin- 
down luminosity after r, Lj oc f" 2 , might agree with that of 
the prompt non-thermal emission flux, F oc f" 2 °, shown in 
Figure I2 Assuming / ~ 10~ 2 , we obtain B ~ 10 16 G and 
Pq r-j 10 ms in order that Lj ~ 10 45 erg s -1 and r ~ 10 3 s are 
reproduced. In this case, the rotational energy Ifl 2 /2 would 
be overwhelmed by the magnetic energy (B 2 /87r)(47r^ s 3 /3) 
at t ~ 10 4 s. This timescale is comparable to the transition 
time of the decay index of the prompt non-thermal emission 
from F oc f" 2 to Fx oc f . If this scenario was true, the ro- 
tation period would evolve as P ~ 10(f/r) 1//2 ms and become 
P ~ 1 s at t ~ 1 yr. If a giant flare occurs at this magnetar like 
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soft gamma-ray repeater (SGR) 1806-20 dHurlev et al.ll2005l) . 
this would be a clear evidence for our proposal. 

The emission mech anism of the prompt no n-thermal emis- 
sion remains unclear. iGhisellini et al.l (12006b) have drawn the 
de-absorbed vF v spectrum from the Swift Ultra Violet/Optical 
Telescope (UVOT) data and argued that the extrapolation of 
the non-thermal spectrum with a constant low-energy pho- 
ton index to the optical band joins with the detected ther- 
mal optical flux. This implies that the non-thermal spec- 
trum does not extend below ~ 10 15 Hz. If this cutoff is 
due to the synchrotron self-absorption in the emitting shell 
with the isotropic kinetic energy L ISQ ~ 10 47 erg s" 1 and the 
bulk Lorentz factor To ~ 5, the emission radius of the non- 
thermal emission is estimated as ro ~ 10 12 cm. It is less 
than the deceleration radius of the afterglow shell estimated 
by 7"dec — 4ro 2 c?d ec ~ 2 x 10 16 cm in our jet model. This im- 
plies that the prompt non-thermal e mission may be produced 
by the internal dissipation of the jet. Ghisellini etldl (l2006bl) 
have shown that a synchrotron inverse-Compton model could 
reproduce the prompt non-thermal emission, the anomalous 
X-ray afterglow, and the UV/optical thermal component with 
using parameters similar to those we suggest in this paper, 
r = 5, 6*o = 0.2, and r = 7 x 10 11 cm. 

The origins of the thermal components of this event are 
also unclear. A possible candidate i s the emission from 
the expanding cocoon a s suggeste d bvlFan et al.l (1200 6) and 
IGhisellini etafl (l2006bl) (see also iRamirez-Ruiz et al.ll2002l) . 



Within the hot jet model discussed in § [4] the energy de- 
posited into the cocoon by the breakout time is estimated 
by Ljtbi ~ 3 x 10 47 erg, which is smaller than the energy of 
the detected thermal emission > 10 49 erg. The rest energy 
in the cocoon is estimated by Mc 2 9 hz 2 /A ~ 10 49 erg, so that 
the cocoon expansion velocity becomes j3 c ~ 10 _1 outside the 
star. The cocoon shell catches up with the afterglow shell at 
t ~ '"dec/c/?e ~ 10 7 s, and thus the cocoon energy will not af- 
fect the radio afterglow. The thermal components are detected 
in the X-ray band at t < 10 4 s and in the UV/optical band at 
10 4 < t < 10 5 s. To explain such a behavior is an interesting 
future problem. 
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TABLE 1 

Spectral parameters of the non-thermal component of GRB 060218 (1) 



f(s) time width (s) a B f3 B E p (keV) \ 2 dof Fis-i50kcV (10 12 ergcm 2 s ') 



219.5 81.5 " L25 -0Q5 " 3 - 8 -62 34 -l9 232J 292 5170 ^40 

332.0 30.0 — 1.47+gg| -8.4' 23+J 156.7 190 6650+'™ 

449.0 45.0 " 1 - 08 -0 06 ~ 2A *%2 H -3 9 239 A 297 9370™'" 

554.0 60.0 " L05 -oI "2^| 10!' 287.8 385 5400+f™ 

674.0 60.0 -1.07X 13+! 391.4 403 4870+ 



-0.06 -0.6 J -4 

,,+0.05 t a+0.3 n+3 



o+^ : 

J14.0 60.0 -1.07+""5 -2.5+Sj 6.5+11 388.5 429 3450+'' 



794.0 60.0 -1-04™"= -2.6+IM 93 390.7 409 3400; 



j_ 15 JOO.J ~>^y jtju_22Q 

1114.0 140.0 -1.17+""* -3.0+"- 2 , 5.5+"^ 629.0 573 1903+J 2 



1404.0 150.0 - 1 - 31 !ooS - 4 -°^2 3 ' 7 -0? 576 - 8 546 89 °ijl 

1704.5 150.5 -l.l3f -2.6™ :2 2.6+j> : ! 438.9 486 830+ 

2005.0 150.0 

2305.0 150.0 



-1.1+J! -2.6+| : | 2.1+f : | 408.6 437 590+{|§ 
-1.1+°| -2.53^-^ 1.3+K2 375.8 393 563 +72 

2593.5 138.5 -l.O+^l -2.8^'j 1.6+° : > 300.9 344 303™ 
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NOTE. — These values are calculated from the joint fit of the BAT and XRT data by the absorbed Band function plus the 
blackbody function. Errors are for 90% confidence. 



TABLE 2 

Spectral parameters of the non-thermal component 
OF GRB 060218 (2) 



t (s) time width (s) a c E p (keV) \ 2 dof 



219.5 

332.0 

449.0 

554.0 

674.0 

794.0 

914.0 

1114.0 

1404.0 

1704.5 

2005.0 

2305.0 

2593.5 



81.5 
30.0 
45.0 
60.0 
60.0 
60.0 
60.0 
140.0 
150.0 
150.5 
150.0 
150.0 
138.5 



-1 36+"- 04 
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-0.08 
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'-0.05 
1 40+O.O6 



-1.44; 
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-1 



:i? 8 
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u iD:i 
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NOTE. — These values are calculated from the joint fit of the BAT and 
XRT data by the absorbed CPL model plus the blackbody function. Errors 
are for 90% confidence. 
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FIG. 1 . — The low-energy photon index ag and the high-energy photon index /3g of the prompt non-thermal emission of GRB 060218 as a function of time, 
from the joint fit of the BAT and XRT data by the absorbed Band function plus the blackbody function. 



1CT r 




Time since BAT trigger [sec] 

FIG. 2. — The light curve in 15-150 keV band and the temporal variation of the spectral peak energy of the prompt non-thermal emission of GRB 060218, from 
the joint fit of the BAT and XRT data by the absorbed Band function plus the blackbody function. The dotted lines show the best-fit power-law decays excluding 
the first two time bins F(t) oc r 20 and E p (t) <x r L6 . 
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FIG. 3. — Schem atic picture of the optical vF v flux from the external shock in our jet model (solid line), which does not exceed the detected optical flux shown 
in Ghisellini et al. 1 2006b) and Campana et al. (2006) (dot-dashed line). The outflow in our jet model is decelerated at t^ ec followed by the Blandford-McKee 



evolution F op t oct 1 , begins to expand sideways at tj followed by the evolution F opt oc t 2 , and finally shifts at t s into the Sedov-Tailor evolution F c 



opt i 



